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FOREWORD 


This  report  presents  the  description  of  a  computer  code  which  was  devel¬ 
oped  to  perform  detailed  fatigue  crack-growth  analysis  on  a  cycle-by-cycle 
basis.  This  conputer  code  is  a  two-dimensional  crack-growth  analysis  routine. 
An  inproved  load  interaction  model  which  accounts  for  both  th  j  tensile  overload 
retardation  and  compressive  load  acceleration  effects  of  a  spectrum  loading 
has  been  inplemented  in  this  program.  The  development  effort  of  this  computer 
code  was  under  Air  Force  Contract  F33615-77-C-3121,  Project  2401,  "Structural 
Mechanics,"  Task  240101,  "Structural  Integrity  for  Military  Aerospace 
Vehicles,"  Work  Unit  24010120,  entitled  "Improved  Methods  for  Predicting 
Spectrum  Loading  Effects."  This  contract  was  administrated  by  the  Flight 
Dynamics  Laboratory  of  the  Air  Force  Wright  Aeronautical  Laboratories,  Wright- 
Patterson  Air  Force  Base,  Ohio.  R.M.  Engle  (AFWAL/FIBE)  was  the  Air  Force 
project  engineer 

This  development  effort  was  conducted  by  personnel  from  the  Fatigue  and 
Fracture  Mechanics  Group,  Dynamics  Technology,  Structure  Systems,  under  the 
direction  of  George  Fitch,  Jr.,  supervisor,  Joseph  S.  Rosenthal,  manager,  and 
Dr.  Leslie  M.  Lackman,  director.  James  B.  Chang  was  the  program  manager  and 
principal  investigator.  Edward  Klein  was  the  original  developer  of  this  com¬ 
puter  code. 
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Section  I 


INTRODUCTION 


This  document  describes  a  conputer  program  which  was  developed  to  per¬ 
form  detailed  fatigue  crack  growth  analysis  on  a  cycle -by- cycle  basis.  The 
fatigue  crack  growth  prediction  method  incorporated  in  this  conputer  code  is 
an  iitproved  methodology  based  on  state-of-the-art  technologies.  The  develop¬ 
ment  effort  and  correlation  of  the  analytical  predictions  to  the  test  tiata 
were  documented  in  reference  1. 

The  baseline  fatigue  crack  growth  rate  equation  chosen  in  the  program 
was  the  modified  Walker  equation  (reference  2)  with  variable  threshold  stress 
intensity  factor  range.  The  generalized  Willenborg  retardation  model  (refer¬ 
ence  3)  and  the  Chang  acceleration  scheme  (reference  4)  were  selected  to 
handle  the  con?)lex  load  interaction  effects  in  a  spectrum  loading  which 
include  tensile  overload  retardation,  and  faster  crack  growth  caused  by 
the  existence  of  the  coitpressive  load  in  cotrpress  ion- tens  ion  (negative  stress 
ratio)  cycles. 

Reduction  of  the  overload  retardation  effect  when  the  overload  is  fol¬ 
lowed  by  a  compression  load  is  also  accounted  for  using  a  reduction  of  the 
effective  overload  plastic  zone  size  approach  formulated  by  Chang  (reference 
1) .  The  linear  approximation  method  proposed  by  Vroman  which  was  widely  used 
in  existing  computer  programs  such  as  EFFGRO  (reference  5)  and  CRACKS  (refer¬ 
ence  6)  was  adopted  as  the  damage  accumulation  scheme.  To  account  for  the 
shape-change  effect  of  a  part-through  crack  (PTC) ,  such  as  the  surface  flaw 
and  a  single  corner  crack  at  a  fastener  hole,  the  two-dimension  (2-D)  crack 
growfth  analysis  approach  used  by  many  existing  crack  growth  codes,  including 
Rockwell's  FLAGRO  (reference  7)  and  NASA's  CRACK  (reference  8),  was  adopted 
in  this  computer  program.  The  2-D  crack-growth-analysis  approach  assumes  that 
growth  in  the  two  principal  directions  of  a  part-through  crack  is  a  function 
of  the  stress  intensity  factors  at  the  extreme  points  in  each  direction.  Tl;e 
CRKGRO  program  also  provides  the  1-D  crack-growth-analysis  option  for  the  nart- 
through  cracks. 

A  collection  of  stress- intensity- factor  solutions  for  various  through- 
crack  (TC)  and  PTC  configurations  has  been  incorporated  into  this  program 
through  a  CRACK  LIBRARY  module  which  consists  of  10  subroutines,  each  con¬ 
taining  a  specific  stress -intensity- factor  solution  for  a  specific  crack  geom¬ 
etry.  There  are  eight  additional  dummy  routines  stored  in  the  program,  which 
provides  the  user  the  capability  for  adding  new  stress- intensity- factor  solu¬ 
tions  for  the  crack  geometries  needed  to  be  considered. 

The  program  provides  the  option  for  counting  the  cycles  for  spectnmi 
loadings  through  the  range-pair  counting  routine  built  into  the  program.  Die 
program  also  provides  the  option  to  perform  parametric  studies  on  parameters 
which  dominate  the  degree  of  damage  such  as  the  design  limit  stresses. 
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The  QRKGRO  program  also  provides  four  options  for  presenting  the  crack 
growth  history  in  graphical  plotting  format.  These  options  are  (1)  crack  size 
versus  number  of  flights,  a  versus  N,  (2)  crack  growth  rates  versus  number  of 
flights,  da/dF  versus  N,  (3)  crack  growth  rate  versus  crack  sizes,  da/dF  ver¬ 
sus  a,  and  (4)  crack  growth  rate  versus  the  maximum  stress  intensity  factor 
per  flight,  da/dF  versus  Kmay.  For  a  parametric  study,  a  maximum  of  seven 
curves  can  be  plotted  on  one  chart  in  order  to  provide  the  user  a  clear 
picture. 
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Section  II 


TECHNICAL  DISCUSSIO'l 


The  technical  approach  used  in  this  conputer  program  is  applicable  pri¬ 
marily  for  metallic  structures  which  contain  cracks  or  crack- like  flaws  sub¬ 
jected  to  cyclic  loadings.  The  crack  growth  analysis  methodology  was  developed 
based  on  linear  elastic  fracture  mechanics  (LEFM)  principles;  i.e.,  the 
range  of  crack-tip  stress  intensity  factor, AK,  is  the  controlling  parameter 
for  characterizing  the  cyclic  crack  growth  rates. 

This  section  includes  a  discussion  of  each  of  the  following  elements 
which  form  the  overall  analysis  methodology  incorporated  in  the  program; 
fatigue  crack- growth- rate  equation,  load  interaction  model,  damage  accumula¬ 
tion  scheme,  and  stress- intensity-factor  solutions. 


FATIGUE  CRACK-GR(MH-RATE  EQUATION 

For  constant- anplitude  loadings,  the  baseline  fatigue  crack- growth- rate 
equations  used  in  this  program  are  the  modified  Walker  equation  (reference  2) 
for  positive  stress  ratios  and  the  Oiang  equation  (reference  4)  for  negative 
stress  ratios.  In  mathematical  forms,  they  can  be  expressed  as  follows: 


For  4K  >  4K^,  R  ^  0 

da/dN  =  C  [4K/(1-R)^''"]^, 


R  4  r""  ^ , 
cut 

R  >  R""  „ , 
cut 


For  AK  >  AK  ,,  R  <  0 
th 


da/dN  =  C  [(1  *  R  K^]",  R  ^  R'ut  ’  ^  ^ 


r<r  ,.,r  =  r  „ 

cut  cut 


For  AK  4  AK. 

tn 

da/dN  =  0 

where  C  and  n  are  the  growth  rate  constants;  m  is  the  stress-ratio  collapsing 
factor,  q  is  the  acceleration  exponent,  and  R*^^  are  the  cutoff  values  for  the 
stress  ratios,  either  positive  or  negative. 
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The  threshold  stress- intensity- factor  range  is  determined  by 

-  (1  -  W  '‘Ky, 

o 

where  threshold  value  of  the  stress- intensity- factor  range  obtained 

from  R  ■  (T constant  anplitude  tests;  A  is  an  empirical  constant  determined  from 
constant-anplitiide  test  data  with  various  stress  ratios. 

A  detailed  description  on  the  procedure  for  the  determination  of  the 
crack-growth-rate  constants  is  in  Appendix  A. 


LOAD  INTERACTICN  MODEL 

Vatioiis  load  interaction  effects  on  the  crack  growth  behavior  under  spec¬ 
trum  loadings  have  been  observed.  The  important  effects  are: 

1.  Tensile  overloads  cause  significant  retardation  of  the  crack  growth. 

2.  Conpressive  loads  accelerate  crack  growth  rates.  Furthermore,  a 
compressive  load  uimediately  following  the  tensile  overload 
reduces  the  retardation  effect  introduced  by  the  tensile  overload. 


TENSILE  OVERLOAD  RETARDATION  MODEL 


To  account  for  the  retardation  effect,  the  generalized  Willenborg  model 
(reference  3)  is  adopted  in  this  program.  The  generalized  Willenborg  model 
can  be  written  in  the  following  form: 


(K  )  =  K  -  (l 

max  eff  *>  max  V 

max  L  ' 

(K  .  )  =  K  -  ( 1 

^  nun  eff  «  .  max  \ 
min  L  N 


da \l/2 
^ol) 


♦  = 


1) 


where  is  the  stress-intensity-factor  corresponding  to  the  maximum 

remotely  applied  stress,  is  the  stress- intensity- factor  corresponding  to 
the  maximum  stress  of  the  overload,  da  is  the  incremental  growth  following 
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the  overload,  ZqL  is  the  overload  interaction  zone  size.  Rgo  is  the  overload 
shutoff  ratio  which  is  defined  as: 


/K, 

^  max/  ' 


max 


For  spectrum  loading,  the  effective  stress- intensity- factor  range  and 
effective  stress  ratio  which  are  used  in  OIKGRO,  are  expressed  in  terms  of  the 
maximum  and  minimum  effective  stress  intensity  factors  as  follows: 

eff  max  eff  min  eff 

R  =  (K  .  )  ^7(K  ) 

eff  min  eff  max  eff 

In  the  load-interact ion-account ed-f or  option,  the  program  utilizes  the  follow¬ 
ing  equation  to  account  for  tensile  overload  retardation  effect: 


For4K  >4K^,  ,  R  >  0 
eff  th  eff 


R 

cut’ 

•^eff 

eff 

+ 

+ 

\ut’ 

\ff  " 

R  , 
cut 

For  4K  j-r  <  -4k  , 
eff  th 


da/dN  =  0 

where  C,  n,  m  and  are  the  same  crack- growth -rate  parameters  described 
under  "Fatigue  Crack -Growth -Rate  Equation."  The  threshold  values  of  the 
stress -intensity- factor  range  are  also  identical  to  those  used  in  the  constant - 
amplitude  cases. 


COMPRESSIVE  LOAD  ACCELERATICW  MODEL 

If  the  effective  stress  ratio  is  negative;  i.e.  ,  Rg££  <  0,  the  Ch;ing 
negative  stress  ratio  equation  is  used  in  this  program,  which  accounts  for 
the  compressive  load  acceleration  effect: 

da/dN  =  C[(l  +  (K  )  R  '  ^cut’  \ff  "  ‘^eff 

eff  max  eff  eff 


eff  cut  etf  cut 


where  q  is  the  acceleration  index  determined  from  test  data  generated  for  a 
specific  negative  stress  ratio  (R  <  0)  and  its  R  =  0  counterpart. 
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The  reduction  of  the  overload  retardation  effect  caused  by  a  compressive 
spike  load  inmediately  following  the  tensile  overload  is  accounted  for  by 
OIKGRO  through  an  effective  overload  interaction  zone  concept  proposed  by 
Chang  (reference  1) .  The  effective  overload  interactive  zone  is  defined  in 
terms  of  the  negative  effective  stress  ratio  (Rg££  <  0)  as; 


(Z  T  )  *  (1  R  T  )  »  R  ££  —  R  ^  »  R  rr  “  R  rr 

oL  eff  ^  eff  oL-^  ’  eff  cut'  eff  eff 


R  cc  ^  R  4.  >  R  cc  ”  R  ^ 

eff  cut  eff  cut 


where  Zq^  is  the  plastic  zone  size  introduced  by  the  tensile  overload. 


In  CRKGRO,  the  plane  strain  plastic  zone  size  is  used  if  the  stress  inten¬ 
sity  factor  at  the  maximum  depth  for  a  part-through  crack  is  to  be  calculated. 
The  plane  stress  plastic  zone  size  is  used  at  the  length  direction  for  TC's 
and  PTC's.  The  plane  stress  and  plane  strain  plastic  zone  sizes  are: 


(Z  1 

oL  plane  strain 


CZ  ■) 

^  oL  plane  stress 


67r 


lit 


2 


where 


is  the  material  tensile  yield  strength. 


DAMAGE  ACCUMULATION  SCHEME 


The  Vroman  linear  approximation  method  has  been  incorporated  into  this 
computer  program  as  the  damage  accumulation  scheme.  The  following  paragraphs 
briefly  describe  the  method. 

For  a  given  load  spectrum  as  shown  in  table  1,  the  Vroman  damage  accumu¬ 
lation  scheme  proceeds  by  considering  a  load  step  (i)  and  using  and 

<Tjnin£  calculate  (da/dN)^.  The  value  of  (0.01a)/(da/dN) ^  is  then  compared 
to  N£,  where  "a"  is  the  instantaneous  crack  size.  If  (0. 01a)/(da/dN) i  is 
greater  than  then  the  crack  growth  for  that  particular  load  step  is 
A  a  =  N£  x(da/dl'l)£,  the  crack  has  then  grown  from  "a"  to  (a  +Aa),  and  the 
program  proceeds  to  the  next  load  step. 

If  (0.01a)/Cda/dN)£  is  less  than  or  equal  to  N^,  the  crack  size  will  be 
(a  +  0.01a),  and  this  load  step  is  reexamined.  This  process  continues  with 
(0.01a)/(da/dN)i  being  conpared  to  the  remaining  cycles  in  the  step.  When  all 
load  steps  in  the  block  or  flight  have  been  examined,  the  program  then  proceeds 
to  the  first  step  of  the  next  block  (or  flight)  and  continues. 
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TABLE  1.  A  TYPICAL  STRESS  SPECTRUM  t'ABLE  fSCHBl^TIC) 


Max  stress 

Min  stress 

No.  of  eye/ 
block  (.flight) 

a  * 

min^ 

or 

max. 

‘^min. 

"2 

** 

°^max. 

X, 

3 

• 

• 

• 

®^max. 

‘^min. 

1 

1 

N. 

1 

STRESS- INTENSITY- FACTOR  SOLUTIONS 

Two  major  types  of  crack  geometry  are  considered  in  this  program;  TC  and 
PTC.  Crack-tip  stress -intensity- factor  solutions  for  commonly  detected  TC 
and  PTC  are  incorporated  into  the  CRACK  LIBRARY  module  of  this  program.  The 
CRACK  LIBRARY  module  consists  of  separate  subroutines,  each  containing  one  set 
of  stress- intensity- factor  solution. 

For  PTC’s,  this  program  accounts  for  the  shape-change  effects  through 
the  2-D  crack-growth-analysis  approach.  The  2-D  crack-growth-analysis  approach 
assumes  that  growth  in  the  two  principal  directions  of  a  PTC  can  be  character¬ 
ized  by  the  stress- intensity-factors  at  the  extreme  points  of  each  direction. 

In  constant -anplitude  loadings  for  exanple,  the  crack  growth  rates  at  the  two 
extreme  points,  A  and  B,  of  a  surface  crack  as  shown  in  Figure  1  are  calculated 
by  CRKGRO  using  the  following  set  of  eqiiations: 

l-m^  n 

da/dN  »  C^  [4K^/(1  -  R)  ] 

1 

dc/dN  =  Cg  [4Kg/(l  -  R) 

where  C^,  n^,  m^  and  Cg,  n_,  nig  are  the  material's  crack- growth- rate  param¬ 
eters  along  the^depth  and  the  length  directions,  respectively*  R  is  the  cyclic 
stress  ratio,  and  andAKg  are  the  stress- intensity- factor  range  at  the 
maximum  depth  and  length  points,  respectively. 

The  stress  intensity  factors  at  the  maximum  depth  point  A,  and  the  maxi¬ 
mum  length  point  B  built  into  CRKGRO  are  prepared  using  the  compound  solution 
format.  In  general,  these  can  be  expressed  as:  . 


where  a  is  the  depth,  c  is  the  half-length  for  a  surface  crack,  t  is  the  thick¬ 
ness  of  the  structure,  b  is  the  half -width  of  the  structure,  Q  is  the  shape 
factor,  and  F;^  and  Fg  are  the  geometrical  magnification  factors  for  the  maxi¬ 
mum  depth  point  A  and  the  maximum  length  point  B  as  sho\m  in  Figure  1. 

In  CRKGRO,  the  gecanetrical  magnification  factors  are  in  a  polynomial  for¬ 
mat.  For  the  surface  crack  shown  in  figure  1,  the  geometrical  magnification 
factors  are  derived  from  the  general  surface  crack  stress- intensity- factor 
solution  proposed  by  Newman  (reference  9).  At  point  B,  the  solution  was 
derived  from  0  =  10“  as  suggested  by  Newman  (Reference  12) . 


8 


9 


At  point  A  (0  =  90®] 


Fa  '  [lo:2°-:'a-/c)'  -  “•=‘‘](f)' 


[< 


+  0.5  - 


(0.65  +  a/c) 


14  (1  -  a/c)^'^]  (|) 


At  point  B,  (0  =  10“) 


Fg  =  ^1.13  -  0.09(a/c)  + 


+  jo. 


5  - 


(0.65  +  a/c) 


+  14  (1  -  a/cl 
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(f) 


! 


1.1  +  0.35(a/t)' 


The  shape  factor  Q  used  in  QUCGRO  is  also  in  a  closed- form  solution  format 
vdiich  was  formulated  by  Newman  (reference  10)  as: 

Q  =  (1  +  1.464  (a/c)^'^^) 


A  collection  of  stress -intensity- factor  solutions  for  various  TC  and  PTC 
configurations  has  been  incorporated  into  CRKGRO  through  a  CRACK  LIBRAF.Y  mod¬ 
ule  which  consists  of  separate  subroutines,  each  containing  one  stress- 
intensity-  factor  solution.  A  crack  code  system  is  used  in  CRKGRO.  ^igure  2 
shows  the  crack  geometries  with  the  corresponding  assigned  code  number.  For 
exanple,  crack  code  1010  is  the  surface  flaw  and  2010  is  the  centered  TC. 

Each  PTC  subroutine  has  two  sets  of  solutions,  K^  and  Kg.  The  stress  intensit)’ 
factor  for  a  shallow  crack  (a/c  <  1)  and  for  a  deep  crack  (a/c  >  1)  are  inclu¬ 
ded  in  the  same  subroutine.  Ten  stress -intensity- factor  solutions  have  been 
incorporated  in  the  CRACK  LIHIARY  module.  Appendix  B  contains  these  solutions. 


For  part- through  cracks,  the  CRKGRD  program  also  provides  the  one¬ 
dimensional  (1-D)  crack  growth  analysis  options.  For  the  1-D  option,  only  the 
stress  intensity  factor  at  the  maximum  depth  of  a  PTC  is  calculated.  The  aspect 
ratio  (a/ 25-)  is  assimed  to  be  constant  through  the  whole  period  of  the  growth 
of  the  PTC,  i.e.,  the  shape  change  of  a  PTC  is  not  accounted  for  in  the  crack 
growth  analysis. 


The  1-D  stress  intensity  factO’*  solutions  for  varioiis  part- through  cracks 
such  as  the  surface  crack,  the  edge  comer  crack,  -are  presented  also  in 
Appendix  B.  These  solutions  are  somewhat  different  than  its  2-D  option  coun¬ 
terpart  due  to  the  fact  that  the  shape  change  effect  is  not  accounted  for  in 
the  1-D  option,  they  were  formulated  using  the  compound  solution  technipue  with 
the  average  value  using  for  geometry  correction  factor. 
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CODE  NO. 


DESCRIPTION 


GEOMElTtV 


SURFACE  CRACK.  CENTERED 


ONE  CORNER  CRACK  FROM 
CENTERED  OPEN  HOLE 


TWO  CORNER  CRACKS  FROM 
CENTERED  OPEN  HOLE 


ONE  CORNER  EDGE  CRACK 


THROUGH-CRACK,  CENTERED 


ONE  THROUGH-CRACK  FROM 
CENTERED  OPEN  HOU 


TWO  THROUGH-CRACKS  FROM 
CENTERED  OPEN  HOLE 


ONE  THROUGH-EDGE  CRACK 


TWO  THROUGH-EDGE  CRACKS 


ASTM  COMPACT  SPECIMEN 
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Section  III 


PROGRAM  OUTLINES 


This  computer  code  is  called  a  detail  crack  growth  analysis  program  and 
its  identification  is  CRKGRO.  A  program  flowchart  is  shown  in  figure  3. 

CRKGRO  consists  of  34  subroutines,  of  which  18  are  basic  to  the  stress- 
intensity-  factor  calculations  identified  as  the  CRACK  LIBRARY.  Currently,  10 
stress-intensity-factor  solutions  have  been  incorporated  into  CRKGRO. 

Appendix  B  shows  these  10  stress-intensity-factor  solutions.  The  other  eight 
subroutines  are  the  dummy  subroutines,  in  which  a  new  stress-intensity-factor 
solution  can  be  coded  into  the  program. 

1.  dlKGRO  -  overall  supervisory  routine 

2.  CENTER  -  centers  titles  with  assigned  fields 

3.  CICON  -  converts  characters  to  binary 

4.  CMBCD  -  moves  BCD  charac:ers  from  one  FORTRAN  array  to  another 

5.  COCEN  -  converts  binary  to  characters 

6.  CRIT  -  conputes  the  critical  crack  length  without  growing  the 

crack 

7.  CYCCNT  -  performs  range-pair  counting  of  the  spectrum 

8.  GROW  -  supervisory  crack  growth  routine 

9.  GETDAT  •  gets  control  card  information 

10.  INPUT  -  reads  the  input  data,  consisting  of  the  crack- growth -rate 

equation  and  load  interaction  model  parameters,  material 
and  gecmietric  constants,  spectrum  input,  and  output  controls 

11.  LIBBD  -  provides  titles  for  crack  library;  i.e.,  stress- intensity- 

factor  solutions 

1?.  NEWLFN  -  changes  file  name  in  FIT  for  versatile  file  number  input 
when  spectra  are  stored. 

13.  NUMBER  -  converts  a  character  to  an  integer 

14.  OUTPUT  -  prints  the  echo  of  the  input  data  and  control  printout 

options 
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15.  PLOT 


t 


-  plots  the  grid,  labeling,  and  data  points 
PTPARM 


16.  PTPARM  -  directs  plotting  of  crack  growth  analysis 


for  subroutine 


results 


Section  IV 

I 

‘  CRKGRO  INPUT  DATA  DECK 


The  input  data  deck  for  CRKGRO  is  described  in  this  section.  A  brief 
description  of  each  type  of  input  data  card  is  presented  in  the  following  list, 
the  overall  deck  setup  is  shown  in  figure  4,  and  a  detailed  description  of 
each  card  follows 


Card  Description 

la  "TITLE”  keyword 

lb  Title  cards 

Ic  "END"  keyword 

2a  "MATERIAL"  keyword 

2b  Material  description  card 

2c  Crack-growth-rate  equation  constants 

2d  Material  properties 

5a  "THRESHOLD"  keyword 

3b  Delta  K-threshold  and  variable  delta  K-threshold  constant 

4a  "LIMITS"  keyword 

4b  Initial  and  final  crack  sizes  and  stress  ratio  cutoff  values 

5a  "ANALYSIS"  keyword 

5b  Load  interaction  specifications 

5c  Crack  code  and  crack  geometry 

5d  Limit  stresses 

ba  "SPECTRUM"  keyword 

bb  Spectrum  title  card 

6c  Spectrum  scale  factor,  range-pair  counting  option,  and  file 

number 

bd  Keyword  for  spectrum  type,  "MAX -MIN,"  "R-DELTA," 

"MEAN-ALT,"  and  flight  mission  segment  title 
be  Stress  spectrum  and  occurrences 

6f  "END"  keyword 

6g  "END  SPECTRUM"  keyword 

6h  Mission  mix 

6i  Number  of  flights  per  block 
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Card 


7a  "OUTPirr*  keyword 

7b  Print  and  plot  options 

7c  Plotter  type 

7d  Plot  types  and  scaling 

8  "END  DATA"  keyword 
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CARO  8 

END  DATA  CARO 


CARDS  7A-7D 
OUTPUT  CAROS 


Figure  4.  CRKGRO  Input  Deck  Setup 


INPUT  DATA  CARDS  NO.  lb 


Description:  Title  cards 
Format  and  Example: 


Column  1  72  80 


TITLE(l)  through  TITLE  (18) 

BENCHMARK  443  FOR  A  DOUBLE  CORNER  CRACK 

Field  Contents 


TITLE  (I), 

I  =  1,  18 


Any  alphanumeric  information  which  the  user  desires  to 
input  for  problem  identification 


Remarks :  There  is  no  limit  to  the  number  of  title  cards  input. 

All  title  cards  are  printed  on  the  first  page  of  output, 
but  only  the  first  title  card  is  printed  on  succeeding 
pages  and  all  plots. 
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INPUT  DATA  CARD  Ic 


Description:  Termination  of  title  card  input 

Format  and  Example: 

Column  1  34 


OPTION 

END 

Field  Contents 

OPTICW  "END"  keyword,  beginning  in  column  1  for  terminating  title 

card  input 
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INPUT  DATA  CARD  NO.  2a 


Field  Contents 

OPTICW  "MATERIAL"  keyword,  beginning  in  column  1. 

Remarks :  This  card  initiates  the  input  of  the  material  description, 

crack-growth-rate  constants,  and  static  material  properties. 
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INPUT  DATA  CARD  NO.  2b 


Description:  Material  description 
Format  and  Example: 

Column  1  60  70  80 


MATID(l)  through  MATID(7) 

BISLP 

2219-T851  ALUMINUM 

Field 


Contents 


MATID(I) , 

I  =  1,  6 

BISLP 


Any  Alphanumeric  information  which  the  user  desires  to 
input  for  material  identification.  60  columns  can  be 
used  for  input 

Option  to  use  single  or  bi-slope  crack  growth  rate 
curves 

Blank  -  analyze  with  single-slope  rate  curve 
"BISLOPE"  -  analyze  with  bi-slope  rate  curve 
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INPUT  DATA  CARD  NO.  2c 


Description:  Crack- gTOVrt.h- rate  constants 

Format  and  Example: 


Column 


Field 

CSUBC 

EXPNC 

CSUBA 

EXPNA 

EXPM 

EXPQ 


1  10  20  50  40  50 


CSUBC 

EXPNC 

CBUBA 

EXP.NA 

EXPM 

z 

5.066  E-10 

3.83 

5.066  E-10 

3.83 

0.60 

7 

60  80 


r 

EXPQ 

L 

1.0 

Contents 

Crack- growth- rate  equation  coefficient  C  used  in  dc/d\ 
equation 

Crack -growth -rate  equation  exponent  n  used  in  dc/dN 
equation 

Crack-growth-rate  equation  coefficient  C  used  in  da/dN 
equation 

Crack -growth -rate  equation  exponent  n  used  in  da/dN 
equation 

Positive  stress  ratio  collapsing  factor  m  used  in  the 
modified  Walker  equation 

Negative  stress  ratio  acceleration  index  q  used  in  the 
Chang  equation 


Remarks :  1.  CSUBC  and  CSUBA  must  be  in  ksi  units. 

2.  If  the  "BISLOPE"  option  is  used  (see  card  2b)  an  extra 
input  card  is  required  for  the  lower  part  of  the  rate 
curve.  This  card  is  similar  to  card  2c,  where  CSUBCL, 
EXPNa,  CSUBAL,  EXPNAL,  EXPML,  EXPQL  is  as  defined 
above  but  for  the  lower  part  i.e.  region  I  of  the  da/dN 
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versus  M  rate  curve,  with  the  additional  input  values 
of  TRANSL  and  ATLEV  in  columns  61-65  and  66-72, 
respectively . 

TRANSL  -  M  at  transition  from  upper  to  lower  curve 
(e.g.:  3.7) 

ATLEV  -  the  level  of  da/dN  at  transition  from  upper  to 
lower  curve  (e.g.:  1.0  -07) 

INPUT  DATA  CARD  2d 


Field  Contents 

CKC  Plane  stress  fracture  toughness  K(-(ksi  nT in.)  used  for  the 

through-crack  instability  criterion 

AKIC  Plane  strain  fracture  toughness  K,p(ksi NTiii. )  used  for  the 

part -through- crack  instability  criterion 

SIGMAY  Material  yield  strength,  or^yCksi) 
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INPUT  DATA  CARD  3a 


Contents 

"THRESHOLD"  keyword,  beginning  in  column  1 

This  card  initiates  the  input  for  the  AK  threshold  option. 


Field 

OPTION 

Remarks: 
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m 


INPUT  DATA  CARD  3b 


Description:  Delta  K- threshold  and  Constant 
Fbrmat  and  Exanple: 

Colmm  1  10  20  80 


DELm 

THA 

2.5 

1.0 

Field  Contents 

DEUCIH  The  fixed  threshold  value  of  AK  at  R  =  o, 

AK  .  (Ksi  ^/ln.) 

THA  The  variable  threshold  constant 


Remarks:  AK^  =  AK^^  (1-AR) 

o 
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INPUT  DATA  CARD  4a 


Descriotion:  Limits  keyword  card 


Field  Contents 

OPTICW  "LIMITS”  keyword,  beginning  in  column  1 

Remarks :  This  card  initiates  the  input  of  the  initial  and  final 

crack  sizes  and  stress  ratio  cutoff  values. 


INWr  DATA  CARD  4b 

Description:  Initial  and  final  crack  sizes  and  stress  ratio  cutoff  values 
Format  and  Example: 


Column 


Field 

CINIT 

CF 

AINIT 

AF 

RCUT 

RCUTN 

Remarks : 


1  10  20  30  40  50  60 


CINIT 

CF 

AINIT 

AF 

RCIH' 

RCUTN 

0.06 

0.90 

0.06 

0.5 

0.75 

-0.50 

Contents 


Initial  crack  length  Cin.) 

Final  crack  length  (in*) 

Initial  crack  depth  (in.) 

Final  crack  depth  (in.) 

The  cutoff  value  of  the  positive  stress  ratio  "R+,"  above 
which  the  material  is  assumed  to  have  no  stress  ratio  layering 
effect  on  the  crack  growth. 

The  cutoff  value  of  the  negative  stress  ratio  "R-,"  below 
idiich  the  material  is  assumed  to  have  no  acceleration  effect 
on  the  crack  grc^h 

1.  If  CF  =  0,  or  blank,  CF  will  be  set  equal  to: 

a.  the  half  width  -  for  surface  and  center- thru 
cracks 

b.  the  width  -  for  edge  cracks 

c.  the  (half  width-radius)  -  for  cracks  at  holes. 

Analysis  will  terminate  when  either  the  critical 
crack  size  occurs  or  the  crack  length  is  equal  to  CF. 

2.  If  AINIT  =  0.  or  blank,  and  the  crack  type  is  a  part 
through  crack,  the  execution  will  terminate. 

3.  If  AF  =  0.  or  blank,  AF  will  be  set  equal  to  the 
thickness.  Analysis  will  terminate  when  either  the 
critical  crack  size  occurs  or  the  crack  depth  is  equal 
to  AF. 
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4.  RCUT  must  be  in  the  range:  0.^  RCUT-«i  1. 

5.  RCUTN  must  be  in  the  range:  -1.0  RCUTN  i  0. 

INPUT  DATA  CARD  5a 


Field  Contents 

OPTICW  "ANALYSIS"  keyword,  beginning  in  column  1 

Remarks :  This  card  initiates  the  input  of  load  interaction  data,  crack 

code,  crack  geometry,  and  limit  stresses. 


INPin'  DATA  CARD  5b 


Field  Contents 

RET  "LOAD  INTERACTIW  keyword,  beginning  in  column  1 

RETTYP  Load  interaction  application.  If  "NO"  is  input,  no  load 

interaction  will  be  considered  in  the  analysis,  and  the  com¬ 
pressive  minimum  stresses  will  be  set  to  zero.  If  RETTYP  is 
blank  or  "YES,"  overload  retardation  and  ccmpression  acceler¬ 
ation  are  considered.  If  "BOTH"  is  input,  then  both  analyses 
will  be  performed. 

RETDA  Retardation  shutoff  ratio 

Remarks :  The  value  of  RETTYP  must  begin  in  column  21. 


INPUT  DATA  CARD  5c 


Field  Contents 

CODE  Crack  code  consisting  of  four  integers,  beginning  in 

column  1 


DIM  Blank  -  Two-dimensional  analysis 

"ONED"  -  One-dimensional  analysis,  beginning  in  column  5 

OTDTH  Width  of  structure,  W  =  2b  (in.)  for  non-edge  cracks 

W  =  b  (in.)  for  edge  cracks 

T  Thickness  of  structure,  t  (in.) 

RADIUS  Radius  of  open  hole,  r  (in.) 

NBRK  Control  for  transition  from  a  part- through  crack  to  a 

through  crack 

For  NBRK  =  0,  breakthrough  occurs  when  a/t  =  1.  For 
NBRK  =  1,  breakthrough  occurs  at  a/t  =  0.5  [1/0.86  (a/lc)] 

Remarks :  1.  Refer  to  figure  2  for  a  description  of  crack  codes 

available  in  the  current  version  of  CRKGRO. 


INPUT  DATA  CARD  5d 


Description:  Limit  stresses  and  analysis  control  parameters 

Format  and  Example: 

Column  1  10  20  30 


NLIM 

SIGLIMCD 

INSTAB 

1 

3 

10. 

BY  MAXIMUM 

1  10  20  30  40 


SIGLIM(2) 

SIGLIM(3) 

SIGLIMCNLIM) 

.85 

1.2 

Field  Contents 

NLIM  Number  of  limit  stresses,  integer,  right-adjusted 

SIGLIM(l)  Limit  stress  in  ksi 


INSTAB  Control  for  determining  when  instability  is  reached: 

(Starting  in  column  21.) 


"BY  MAXIMUM"  -  by  stress -intensity -factor  value  at 
maxiimm  spectrum  stress 

"Blank"  -  by  stress -intensity- factor  value  at  limit  stress 

SIGLIM(l) ,  Ratio  of  the  ith  limit  stress  to  the  first  limit  stress 
I  =  2,7 


Remarks :  1.  If  NLIM  is  greater  than  1,  then  a  parametric  study  is 

performed  with  the  change  in  limit  stress.  The  spectrum 
stresses  also  change  by  the  ratio  of  limit  stress  change. 

2.  Limit  stress  ratios,  SIGLIM(2-7),  are  input  in  fields 
of  10  on  the  next  data  card.  A  maximum  of  six  limit 
stress  ratios  may  be  input  on  this  card. 


INPUT  DAT/  CARD  6a 


Contents 

"SPECIRIK'  keyword,  beginning  in  column  1 

This  card  initiates  the  input  o£  the  spectrum  data. 


Field 
OPTION 
Remarks : 
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INPUT  DATA  CARD  6b 


Field  Contents 

STITLECI),  Any  alphanumeric  information  which  the  user  desires  to  input 
I  =  1,7  for  overall  spectrum  identification 

Remarks :  Only  information  in  columns  1-70  is  printed. 
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INPUT  DATA  CARD  6c 


Description:  Spectrum  scale  factor,  range-pair  counting  option  and  file 

number  for  stored  spectrum  fif  any) 

Format  and  Example: 

Column  1  10  15  20 


FAC 

NRP 

NFILE 

1.0 

1 

12 

Field 

FAC 

NRP 


NFILE 


Remarks ; 


Contents 

Factor  used  to  scale  the  spectrum 

Control  for  range-pairing  each  spectrum  segment: 

0  -  no  range-pair  counting  operation 
1  -  range-pair  counting  option 

Chit  number  of  file  where  spectrum  is  stored: 

Blank  -  spectrum  will  be  read  in  from  cards 
10-98  -  spectrum  will  be  read  in  from  given  file 

1.  All  maximum  and  minimum  stresses  are  scaled  by  the 
value  of  FAC. 

2.  The  values  of  NRP  and  NFILE  are  integers,  right-adjusted. 
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INPUT  DATA  CARD  6d 


Description :  Spectrum- type  keyword  and  segment  title 

Format  and  Example: 

Column  1  10  70  80 


SPCTYP 

SEGTTL(l)  through  SEGTTLC6) 

MAX-MIN 

AIR-TO-GROUND  SPECTRUM 

Field  Contents 

SPCTYP  Spectrum-type  keyword  defining  the  following  data  cards  as 

one  of  the  following: 


"MAX-MIN”  -  maximim  and  minimum  stresses  are  input 
"R-DELTA"  -  stress  ratio  and  delta  stress  are  input 
"MEAN"  -  mean  and  alternating  stresses  are  input 

SEGTTLCI)  Mission  segment  description 
I  =  1,  6 


Remarks : 


1.  The  keyword  for  SPCIYP  begins  in  column  1. 

2.  The  mission  segment  description  can  be  input  in 
columns  11-70. 

3.  Cards  6d  to  6f  must  be  repeated  for  as  many  segments  as 
there  are  in  the  spectrum  to  a  maximum  of  20  segments. 
There  is  no  limit  to  the  number  of  steps  per  segment, 
although  a  limit  of  3,000  steps  is  established  for  the 
entire  spectrum. 

4.  This  card  should  always  be  a  part  of  the  input  deck 
Cnot  a  part  of  a  stored  spectrm) . 
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INPUT  DATA  CARD  6e 


Description:  Stress  spectrum 

Fonnat  and  Example: 

Column  IS  15  25  35 


SI 

S2 

CYaES 

12. 

-2. 

1. 

25. 

10. 

.01 

6.5 

2.1 

25. 

8.4 

5.2 

500. 

Field 


SI 


S2 


CYaES 

Remarks: 


Contents 


Depending  on  the  spectrum  type,  SI  is  one  of  the  following: 

1 .  Maximum  stress 

2 .  Delta  stress 

3 .  Mean  stress 

Depending  on  the  spectrum  type,  S2  is  one  of  the  following: 

1 .  Minimum  stress 

2.  Stress  ratio 

3.  Alternating  stress 

Number  of  occurrences  for  each  type  of  loading 

1.  The  spectrum  is  considered  to  be  in  ksi  units. 

2.  A  value  for  CYCLES  <  1  will  be  applied  every  (CYCLES)  ^ 
times.  If  CYCLES  =  0.1,  the  loading  will  be  applied 
every  10th  time  that  the  flight  segment  is  repeated. 

3.  A  maximum  of  3,000  steps  is  established  for  the 
entire  spectrum. 

4.  If  the  spectrum  is  in  terms  of  %  design  limit  stress, 
the  scale  factor  (card  6c)  will  convert  the  spectrum 
into  the  correct  stress  level. 
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INPUT  DATA  CARD  6f 


Field  Contents 

ENDSEG  "END"  keyword,  beginning  in  column  1  for  terminating  flight 

segment  spectrum  input 

Remarks :  1,  The  succeeding  data  card  is  either  a  flight-segment-type 

card  6d  or  card  6g. 

2.  If  the  spectrum  is  stored  on  a  file,  then  the  keyword 
"END"  must  also  be  stored  on  the  file  after  each  segment 
and/or  at  the  end  of  the  spectrum. 
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INPUT  DATA  CARD  6g 


Field  Contents 

ENDSPC  "END  SPECTRUM"  keyword,  beginning  in  column  1  for 

terminating  spectrum  input 

Remarks :  This  card  should  always  be  a  part  of  the  input  deck  (not 

part  of  a  stored  spectrum). 
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INPUT  DATA  CARD  6h 


Field  Contents 

MISSION  "MISSION  ="  keyvrord,  starting  in  column  1 

NBLKS  NuiTiber  of  times  the  complete  mission  string  is  repeated 

FACTORi  Number  of  times  the  individual  flight  segment  is  repeated 

SECWENTi  Mission  segment  number  preceded  by  the  letter  "S"  for  segment 

Remarks :  1.  The  exanple  illustrates  an  equation  for  mission  mix. 

The  characters  "MISSION  ="  are  required  to  initiate  the 
mission  mix.  The  parentheses  are  required,  and  the  tenns 
between  them  describe  a  conplete  mission  string.  All 
characters  are  free  format. 

2.  The  maximum  number  of  individual  mission  segments  is  20 

CSl . S20) ,  and  the  maximum  number  of  mission  segments 

mix  is  100  C1S1+ . ). 

3.  No  limit  is  est.’blished  for  NBLKS. 

4 .  Other  examples : 

MISSIOM  =  10000 

MISSION  =  500  (ISl  +  1S2  +  1S3  +  1S4  +  1S5  +  1S6 

1S5  +  1S4  1S3  +  1S2  +  ISl  +  1S2  +  1S3  + 

1S4  +  1S5  +  1S6) 

5.  The  last  entry  must  be  a  "+"  sign  for  continuation  of  a 
mission-mix.  The  end  of  the  mix  is  a  closing  parenthesis 


42 


INPUT  DATA  CARD  6i 


Field  Contents 

NFPB  Number  of  flights  per  block,  integer,  right- adjusted 

Remarks :  1.  This  value  is  used  for  plotting  only. 


INPUT  DATA  CARD  7a 


Contents 

"OUTPUT'  keyword,  beginning  in  column  1 


Field 

OPTION 


Remarks : 


This  card  initiates  the  input  of  printer  and  plotter  controls. 


INPUT  DATA  CARD  7b 


Description: 
Format  and  E 
Column 

Field 

IPRSPC 

NB 

NCRT 

NJDMP 


IPR72 

KPRT 


Print  and  plot  options 


1  5 

10 

15 

20 

25 

30 

35 

40 

45 

IPRSPC 

NB 

NCRT 

NJIM» 

IPR72 

KPRT 

INPRT 

LNBLK 

1 

10 

1 

1 

0 

1 

0 

-25 

Contents 

Control  for  printing  the  spectrum  to  be  used 
1  -  print 

0  -  suppress  print 

Control  for  printing  the  crack  growth  history  in  increments 
of  NB  number  of  blocks.  The  default  value  is  175. 

Control  for  plotting  and  reading  plotting  data: 

1  -  read  plotting  parameters  on  succeeding  cards 
0  -  no  plotting 

Control  for  bypassing  slow-growth  steps: 

^8 

0  -  bypass  steps  with  a  growth  rate  less  than  10  and 

M<.2M 

max 

X 

1  -  retain  slow-growth  steps  in  analysis 

Control  for  limiting  the  line  size: 

1  -  printed  output  will  have  a  72-column  width 
0  -  printed  output  will  have  a  108-column  width 

Control  for  printing  stress  intensity  equations  used 

1  -  print  K  equations 
0  -  do  not  print 
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INPUT  DATA  CARD  7b  (Concluded) 


Field 

INPRT 

LNBLK 


Remarks : 


Contents 


Control  for  printing  growth  history  of  first  block 

1  -  do  not  print 
0  -  print  first  block’s  growth 

Control  for  printing  intermediate  growth  history 

-i  -  growth  of  the  first  ith  steps  of  each  "NB"th  block 
will  be  printed 

+i  -  growth  of  the  last  ith  steps  of  each  "NB"th  block 
will  be  printed 

o  -  all  steps  of  "NB"th  block  will  be  printed 

If  the  spectrum  is  segmented  and  the  first  segment's  number 
of  steps  <LNBIJC,  printing  will  stop  at  the  1st  segment-end 
for  "-i"  or  start  at  the  last  segment-start  for  "+1".  First 
block's  print  will  also  follow  LNBLK 's  option. 
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[NPUT  DATA  CARD  7c 


Field  Contents 

NPLTYP  Plotter  type  for  DISSPLA  processing: 

1  -  SC4020 

2  -  Tektronix 

3  -  Unipost  or  Calconp 

IBAUD  Tektronix  terminal  line  speed  in  characters  per  second 

(used  only  when  NPLTYP  =  2) 

Remarks :  1.  All  values  are  integers,  right -adjusted. 

2.  Card  7c  is  input  only  if  "NCRT"  on  card  7b  is  1. 
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INPUT  DATA  CARD  7d 


Description:  Plot  types  and  scaling  parameters 

Format  and  Example: 

Column 


L  5 

10 

IS 

20 

25 

30 

50 

55 

60 

OPTCD 

0PTC2) 

0PTC3) 

0PTC4) 

SCALE 

(2,1) 

■ 

SCALE 

(1,4) 

SCALE 

(2,4) 

■ 

1 

1 

1 

1 

1 

0 

0 

■ 

0 

1 

Field 


OPTCl-4) 


SCALEC2,4) 


Remarks : 


y  Y  vs  X 


Contents 

Array  specifying  parameters  to  be  plotted: 

1  -  plot 
0  -  no  plot 

OPTCl)  =  crack  size  versus  life  in  flights 

0PTC2)  =  crack  growth  rate  versus  life  in  flights 

0PT(3)  =  crack  growth  rate  versus  crack  size 

0PTC4)  =  crack  growth  rate  versus  maximum  K 

•  1.  -  nisx 

per  flight 

Array  specifying  linear,  semilog,  or  log-log  grid  scaling: 

1.  SCALE (1,*)  =  0  -  X-axis  is  linear 

1  -  X-axis  is  log 

2.  SCALE(2,*)  =  0  -  Y-axis  is  linear 

1  -  Y-axis  is  log 

1.  All  values  are  integers,  right -adjusted. 

2.  For  part-through  cracks,  separate  plots  are  produced  for 
both  crack  length  and  depth. 

3.  Because  one  block  represents  the  whole  segment  series  or 
mission-mix  expression  between  the  parentheses  on  card  6h, 
the  number  of  flights  cn  card  6i  must  reflect  both  the 
multiple  applications  and  segment  repetitions. 


4.  Card  7d  is  input  only  if  "NCRT"  on  card  7b  is  1. 
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INPUT  DATA  CARD  8 


Field  Contents 

OPTION  "END  DATA”  keyword,  beginning  in  column  1 

Remarks :  This  card  terminates  the  reading  of  input  data.  If  no 

input  data  errors  were  encountered,  execution  begins; 

,  otherwise,  the  program  is  terminated. 
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Section  V 


EXAMPLE  CASES 


This  section  presents  two  example  cases  which  are  designed  to  illustrate 
the  capability  of  the  CRKGRO  program  and  exercises  most  of  the  program  options. 
The  first  example  is  the  crack  growth  prediction  for  a  surface  crack  contained 
in  a  finite  width  plate  as  shown  in  figure  1,  subjected  to  spectrum  loadings. 
For  the  purpose  of  illustration,  the  spectrum  selected  for  the  analysis  was  a 
small  block  of  random  flight  spectrum  which  consists  of  57  cycles.  Each 
cycle  has  different  stress  levels.  The  analysis  was  done  by  the  two-dimension 
crack  growth  analysis  option.  Crack  growth  rate  constants  used  in  the  analysis 
were  the  bi-slope  constants.  The  following  data  were  used  in  the  analysis. 

Material: 

2219-T851  Aluminum  Plate 

Region  II  (upper  slope)  crack  growth  rate  constants 

C  =  5,066  X  lO'lO  (in  ksi  unit) 
n  =  3.83 

Bi-slope  transition  point: 

da/dN  =  6  X  10'^  in/cyc  AK  =  5  ksiyiiT 

Region  I  (lower  slope)  crack  growth  rate  constants 

C  =  2.126  X  10'13  (in  ksi  unit) 
n  =  9.23 

Crack  growth  rate  parameters  and  fracture  properties: 


m  = 

0.6 

a  =  1.0 

1) 

o 

2.5  ksi  /liT 

A  =  1.0 

+ 

‘cut 

+0.75 

R'  =  -0.99 
cut 

R  = 
so 

3.0 

O’  =48  ksi 

ty 

Kr  = 
Ic 

45  ksi  yirT 

=  65  ksi  v/TrT 

50 


Plate  dimensions: 


2b  =  6.0  in. 

Initial  crack  sizes: 

a.  =  0.10  in. 

1 


t  =  0.25  in. 


a. /2c.  =  0.5 
1  1 


All  the  inpu"  ecnoes  and  the  nrint-outs  of  the  output  including  the 
graphics  aio  shovni  in  the  next  few  pages.  Brief  descriptions  are  provided 
for  each  page  of  the  output  printouts. 
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they  can  be  included  directly  in  the  crack  growth  analysis 
report , 
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Remarks :  This  is  the  printout  of  the  critical  crack  size  calculation 

based  on  the  design  limit  stress.  The  calculation  is  per¬ 
formed  through  subroutine  CRIT  by  an  iteration  procedure. 
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Remarks :  This  is  the  printout  which  provided  the  detail  data  related 

to  the  cycle-by-cycle  crack  growth  analysis.  Again,  only 
the  first  20  steps  of  the  spectrum  were  printed  out. 
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Remaj'ki;;  This  printout  provided  the  information  on  when  tlie  surface 
crack  became  the  through-the-thickness  crack.  The  total 
cycles  a])pl ied  were  S096  x  57  +  16  x  1  =  290,488  cycles. 
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SAfPLE  CASE  FOR  SURFACE  FLAW 

TFE  CRACK  AT  THE  L'EGINM  IG  0^  tL'>CK  C  =1.711667 
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CRK6R0  PROGRAM  CRACK  GROWTH  ANALYSIS 

SAMPLE  CASE  FOR  SURFACE  FLAW 

LORO  INTERACTION  CONSIDERED 


CRACK  DEPTH  tlNCHES) 


CRKGRO  PROGRAM  CRACK  GROWTH  ANALYSIS 

SAMPLE  CASE  FOR  SURFACE  FLAW 

LORO  INTERACTION  CONSIDERED 


0.0  SOO.O  1000.0  ISC0.02QCn.02500.03000.03SD0.04000.0<SOO.OSOOO.OS500.Q 

LIFE  (FLIGHTS) 
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DC/OF  (INCHES  PER  FLISHU 


CRKGRO  PROGRPiM  CRACK  GROWTH  HNRLYSIS 

SHMPLE  Cf6E  FOR  SURFfiCE  FLRW 

LORD  INTERRCTION  CONSIDERED 


LIFE  (FLIGHTS) 
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CRK6R0  PROGRAM  CRACK  GROWTH  ANALYSIS 

SfiriPLE  CfISE  rOR  SURFACE  FLAW 

LORD  INTERACTION  CONSIDERED 


[iWiI 


CRKGRO  PR06RRM  CRRCK  GROWTH  HNFILYSIS 

SnnPLC  CASE  for  surface  flaw 

LORD  INTERACTION  CONSIDERED 
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The  second  example  is  a  test  data  correlation  case.  Test  data  were 
generated  in  Phase  III  experimc atal  verification  program  (Ref.  12).  The  test 
specimen  in  this  example  case  ;’s  shown  in  figure  5,  which  is  a  2219'T851 
aluminum  center-crack- tens ion  (CCT)  specimen.  It  is  a  standard  ASTM  specimen 
used  for  da/dN  testing.  The  initial  half  crack  length  was  Ci  =  0.148  in, 
which  was  the  measured  dimension  after  the  application  of  20,000  constant 
amplitude  precracking  cycles.  The  test  specimen  was  the  fighter  air-to- 
ground  (A-G)  mission  generated  in  Phase  III  of  this  research  and  development 
work.  Appendix  C  contains  the  peaks  and  valleys  of  the  test  spectrum  in  the 
percentage  of  the  design  limit  stress  (%  «■  format. 


Since  the  crack  was  a  center  through-crack,  crack  code  2010  was  called- 
out.  The  crack- growth -rate  constants  and  load  interaction  model  parameters 
used  in  the  prediction  were  as  follows.  In  this  example,  the  single-slope 
crack  growth  rate  data  was  used. 


th. 


R 

so 


5.066  X  10  (in  ksi  unit) 

q  =  1.0 

3.83 

r""  ^  =  +0.75 
cut 

0.6 

R'  =  -0.75 
cut 

2.5  ksi / in. 

A  =  1.0 

65  ksi  v/in. 

cr  =  48  ksi 
ty 

3.0 

Both  the  load- interaction  and  no  load- interaction  options  were  executed 
in  this  example  in  order  to  illustrate  the  spectnmi  load  interaction  effects . 
For  the  load  interaction  solution  option,  CRKGRO  predicted  the  crack-growth 
life  was  Np  =  19  x  263  +  1,452/19  =  5,073  flights,  or  19  x  4,997  +  1452  = 
96,395  cycles.  Compared  to  the  test  result  (%  =  5,403  flights),  the  predic¬ 
tion  ratio  is  Np/1^  =  0.94.  The  noload  interaction  prediction  was ’Np  ==  3,948 
flights  (75,016  cycles),  compared  to  the  test  result,  Np/N^  =  0.73. 


All  the  input  echoes  and  the  print -outs  of  the  outputs  are  shown  in 
the  next  few  pages. 
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Material 

2219-T85i  aluminum  plate 


Fijmre  5.  Test  Specimen  Configuration 


AJR-TO-GROINO  SPECTRUM*  BASELINE 
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HAXlKLH  KUHBrk  OF  LQtD  BLCCKS  =  200 


AIK-TO-GRCUM  JHrCTRLM*  BASELINE 


UL  iT  id  •-»  C4r^CJf*^»-lfOW^*«<^rSi#^M4r>^#Of**<MrO  f\J  *1  r-<  ir>  ^C-  fO  Al<CgC>jrviC\i<'jK>f^-^K> 

®  r  I  *1 1  r  1 1 3 1 1 1 1 1 1  r  1 1 1 1  r  i  •  i  r  i  r  111*1  1 1  •  1  <11 

^  «iii.  ■T'flO'ifi  vi3a'if)c»r»f0^^CMf^\0r*K)0^^0'U)«-<of0»j0^ifi  ao<N/u;occru®f^coir)'ii^o^r*csi 

<\i  XU.  ooiootaflCi-4’®ir)\oo<M®cOsd'>fluatjvcj\<rir/*-icjoMD^ ff'<r 

o  C5Ui  ••• 

j\  vxy“«^TO K>NO*TNe  xt*-4fox 

(/y  *.4 


o  ^ac(j\  Bs3-x'oooco'X}r^LnK)C'j--rr^®f^-T 

o  cxjoio  inc'C\i«-4«HrOoo^ooocN4<\to<N40CM^e^<4‘vHCa«-4^o<q-iOr-*oo*^r*c>o<V4:>^mroioc3 

o  ac  •••  *•••#••••.•••••••••••••••*••#••••••••••••••• 

fr>  III  I  <=»  f 

«-•  i 

•  X  ooo  oooooooooooo^oooooooooooooocdooooooooooooooo 

<  f^rj ^  r«<MK)h« to XiTJ^tn *-*0"  <^•-llO^~mfOv£rou  r-la:c^J^^ 


o  iDiD  •-4^Groa^^^<4’^ocr\««iinQCoinu3inr^^oac'cg(n^o^o^tof^in-^«^«^^^£^^ 

)  ®®v  ;  p  jo<MCN<Ai) ^^f**^Q»ocjf^cur^K>»^ <r t-4««^f*"C\jsflcr>or^  id 

•-N  X  •#•••♦•♦•••••••••#•♦•♦•••••••♦••••♦••••••••••• 

X  ^ 

X  Kxvjso  K><<r^fOCM\oc>jr»Oso^4r  >iJf^3M4ntf»^.nK)u’)<r  ®<vjp^  O' ocMotocjp-^^ 

^  u  {NiCT  *CM  »0?0vC  ®K>  O'i5<ViK>X'C>C0t^‘X<r  ^'Cg CJinif)®  f^ 

<«)  ^  ••«•••••••••«•••••••••••#••••###••#•»••*•••••• 

c  ^  «■ '®u- cviT  ^^jror-oc'Tuisfl^r* r;cv.<v^r  r  <rK'N/)r*' irt^C4«oc.ino*  r 

^  1  ^  ^  ^  ^  ^ 

X  o 

^  0,  a  >x- 

O  ‘-iJc*  o«jj 

I  •  I  I  I 

'y  ;?  Uiu^iij  uiiJ 

»  l1  t  .— r\  P-. 

^  ^  T>vP  OOfO 

cj  «  \C^ 

^C3  •••#•.. 

>-^  r*  cg^  CNJ O' pH  om  >0 •H IT  P^Oi  X  ^org  ^ i-* OK) IN  pH ro  s3CMu'^  O' ON K)gj 

'Jj 

^  rg  <M<M  cg<'j<MCv/cvi<gcjoj':g<\*cgojfgc>joac\i  cgrgcgcgCMcgcMCVicNirjojcsjcNicgcMcgcgcgCMcvjcMCNjcgrgrg  cjCM 

i  a'*/'a'(7'TJ'r‘crvj'(rcrjN(ro^O'j'(yvcrvo^j'U',riT»'T  tt'C  J'T'0^’T  j'crav0'3^a'C'0'J‘jNcra*'7^ 

^  O  pH  #<»■<»<  iHpi<  pp<  I  I  I  4  >H  ^  ^  pH  pHi^  <— « pp<pH  ^  rHP*** 

r  ••••••••••••«••••••••••••#•••••••••••••••••»•• 


\n  frr*  r^r  >^jr-wr-r-ror^  xf^A-xivc 

04.*w^o  oocioo^crci  o<o<^Ow 

iitii  litiift*  ii»ii 

ijJIaJi^vuuj  ItJ  Ui  I  *-i  ejJ  iaj  tij  \  1  i 

r>i  r  p^■^  a-  r  :^t  .1  ^  r-  ■' 

K5*Hr*p«4X)y-jc\i  O' 

in  <TkoN-<M  <Ts (DpHinm  lO X  ®  f*p-^iCpHAo 


^  c.  \i.N.  X  r  h«  Nr'i;\£  g.!^  a'r 

t**.  o  0  w  c”  Ci  u-#c»  o  w  c^c; »  .  c:  o  wTO 

ll••ll•llll•ll•l  «* 

UjUJlUUUJ(4J‘^'uJUjU^>UU;UJUJ(«J  iu  wJ’iJ 

r^li'  n  l* 

r*-rgp-*a>fOiP  sarvjin  rnr-ro  'x  o  r«H  x 
•.»'XKfl>fOv.3f^pHinot-»»ovON‘ CT'Ln  r^*-< 


H*  (/)  «•••••••••••«••«#«•••••••••»•••••••»•••••••«•• 

^  u 
o  >- 
<  o 

nj  x  c^p-iCj  K)«rto  (T' *3  «Hfg  ro-^  vX)f^  X  O' r>.-^  eg  K' X  O' pH  fvifo  ^  n  XT' 

♦  o  0.  tDiniD  jTu'iu’ikniDuiUi >43'^ 'iJvflvO'fl jnojwoxinxo^  j'<j''T\  j> cr* 

Xi  ui  O' (Tvov  j^'C- J' jNO' vT O' o'dtcr O' O' 0''^<r‘(r  o'jN^xio^o' J'jv  jNtO' O' j't7\o> (J'cyvffs a^lTvu^(^^a^  O' O' J' cT>  J' 

•  ui  ►-  ?^cgej  cjcgejcgcviCMCgo4c\iCSj<>4Cgfg<\«>4<N.rgcgrg<MrgrgrjfgCNjC\icgc><rgcN*cgcNic>jcjcg<\.tMC\j<Ni04'NiC\iCM 

X.  i/j 


86 


STEP  CYCLES  C  CC/uN  DELTCK  KM/X-C  SIGMAX  R  SIGMAX  R 

(EFF)  (EFF) 


r»- C>i  CM  If)  O  O  ■♦•  m  C3  *-«  O  C\i  O  e\J  «  r*  •«•  O  CM  C\J  lO  O  CJ  Cg  CM  <\J  t-l  O  vii  r<  If)  fl- ^ o  o  o  OJ  o  C  J  K ; 

lill  lil  II  lllli  III  I  III  II  II  III  iliii 

(^CMO^>lf)oo«■•-la|f)^OK^•oao•^oo^O^-tf)^\£(Ooo^^•vc^~o•«'(^^f)(^a^  irmino'  * 

oc '4' o  SD  lo  r»  lO  o  o  ifm  <«  «■ «' 00  o  m  r* -4- If)  O' CM  K>  vfl  r<- K)  CO  M3  cc  if> 4- If)  «■  in  JO  4- If)  CM  i> )  <■  o  u '  so 

OOfO^CMsOOCMCViOoaCMCOCMIOsHtoe  OS>-ilomif)s£COh>  ^OmaOsOIO•-•<S4-OaC•HIOI-l4-|^sOroosO<^ 

^CMCMr“lf)in4'CMir»xcO»^<'l4'sO(Ts«-«4’*''‘=Hf)4‘l')4-|nsl)0  ffSsOCMsOci  0'4‘K>  '>CM»-<s£.aor-(Ms-(rt  vC  4- 

^  Oi  ^  ffH  fM 


f04r^<'jxoeoK><TO«-* jojo«^o *fOCNi»o«-^coi-^*-<ojr^ifi  sidk: o^oocf^x-j-vj-ir-M 
fO  a  a.^  O  <V  ^  ^  CJ  ro  C3  ^  CM  O  ^  CV  K)  \0  ^  ^  «H  *0  o  O  ^  O  O  lO  ^  ^  LD  c»  ^  o  CNi  P }  kp  f-t  .Q>  r-l  CsJ 

III  I  I 


ooaoooooooooooooooooooooooooooocsoooooooooooooo 

9^QC««40CL'v-«(T^ro''4>cMNCvo*ocMf^o 

sf)«'00raO^lf)^■^alS•H.HK)«  •a(TS*H«-*lf)^CMO'<^^Os^)r^fMC‘JO«•^<-IOCMOC^•lO■^  vO^-•-•C:^^-K^  J)0  03 

^  •Ht  «4  CM*4ff  II  I  sHCM*^  9  iH 


o4om«H^M34‘io  oo<soscr>>oif)^r»0'so  s'lncMO^msA'^'^gcDsfir-  f^sd  004-0  cM4's£io(Mco(M'osoinin 

I04''fi0(1^4'(l'lf>4‘lf)inf->0'  9'4- ^0^-^-s0Ofs.«4••HOI^t-^s0^-■C0K>O  ®*-«»gin4-  I^CMCOCM®00r>•sO4- 


Il^JWCV«-^rj4‘U#»"“^CS(4' s*J4-s**\s*l—j*)l»3%Sifc^^^l*s*#4- S43r'si-*s*rfi»i  iMr-'s-s#  i 


K)io«oi0  4-o»-iincMsfloocMai  j'4-or>-*-**-<0's0s0f-  ostocMr-Ort 
acM<3if)if)sacMio4'vocr>-iO'«Hjs  000  cm  ®cMif)4'CM*-ieMa"f>4-<jD 

•  «••!••••«••  •<•••••••  •'••••••••• 

ir,c3<Mf"«o<«vK- voiTseo-Kjcr  <rcorjiof<»4-fM»r  9‘c.^a:  tctn 


sf)  r*  co®ac 
o  0000 

I  •  •  I  I 

uj  uijuiu 

Xu-sOlP 
•T  •T'pr-if) 
•H  ^  »ac^ 


« iTh-  r-  sf  X  f*-  r--  IT.  4.  cr®  a'  sx.  sc  r-  M.  r--  f-~ 

O  00  sO  sO  CO*  sJ  ^3^0  COr  C33  c^  c.  I O  C3  C3  sfj  <3>  sO 
UJUJUUJUJ  UULlJtiJUj  luUiUll  ilUJ  IsJUJIiiUJ 'sJUJ 

S'  • ,'  J-  c  ■t^r  fr  <r  IS  •'0»-u'  4-  v  :  >•  — 

in  'f)4-msor^mO'»o<js®in  ®  xcsjm  ®o 
^■if)r*-m»o<Minr»sO®«-<OK)K)K)  ®.-<<T>-<rjcM 


eoo4-sf)h-  «in®r“4-(i''"'»“'ioin>f)'~ 
O  ®  Ln  f"  O  CM  If)  O' I*- ®  CM  r-*  rl  O' fs- 


•  •  •  • 

♦  • 

• 

OCV.C 

cv^ir^ 

? 

8 

CM  IT  4- 

sf.voci'r  1"  rr'-sT 

sT 

si; 

^  X 

L_/0<i0  1^0  oo 

GC 

1  1  •»  1  1  1  1 

1  •  1 

1 

•  1 

*  U  tiJUJ  LJ  UJ  UJ  LwUJ 
T'^  •■*•  .  O'  T  C 

UJUJUJ 

UJ 

UJUJ 

S^  < 

r-®MJP*CM^OK*i 

OCMU' 

X 

•a-  <\J 

J^oXi-^CT) 

sxrOK) 

<■  X 

OCM®  ®r~'»'0®m.HIOr»<'JfO»^»-<®4-i|'0'  CMiMCM4-«-lr^  ff'r'tCM®  »4-s-«t-<lO»-''^in4-CM  => .-I O'M  O' 


CMPfCMCMCMCMCVlCMCMCMCMCMCM  CMCMCMCMCMCMCMCMCMCMCMCM  CMCMCM  CMCM  K)rOK)mK3  lOK)  fOfO  lOlO  K3  rO>0  "1 
O'O'O'  O'CJCMCMCMCMCJCMCJCMCJCMCMCMCMCMCMCMCMCMCMCMCMCiCMCOCMCwC'JCMPJCMC'oCJCMC'xCMCVfMCMCMCMCJ 
4-4-'^-^®  ®  ®if)®®m  -fiui®  ®  ®:nminintf)'n®  ifiT,  uiininuTfiminiD  niDininuTininininn  in 


u^inaOiH ® ‘-C  s®® x  o 
(j\0\(7v  ^ ^vr^g>iys  j'  j'i/%o>(js(yf^ J'ctxtso^j^  (r'^o>C7^  j* 

<M  CM  CM  CM  ^  ^  ^  ^  ^ ^  ^  ^  ^  ^ ^  ^  ^ ^ ^  «r 'T ’T 'T  ^  <9- ^  <9- ^  sT  s^  -T 


87 


AIR-TO-GKCJN  :  S^':CTRUM.  BASELINE 


CCiL, 
00  UJ 


C\)0«H  too  If)  sO  ID 

•  ••  ••••• 

ill  I  I  I  I  f 


o 

h-cc***  inxcsiOiD 

"V 

CM 

s:u. 

<9‘<30ir  mooocMo 

C3 

c  u 

«••  ••••• 

U4  W 

rMlD'T  <\J>OtOoa<M 

VJ 

vH  ^  ^ 

O 

a 

0  0*H  tOfHOOlT) 

o 

OL 

U) 

1 

X 

• 

X 

ooo  0C3000 

r- 

fOCMCr  4}ODXCVj 

CM 

11 

X 

f*-r-oo<MCj 

CM 

o 

Si 

o 

4- O' 00 

</1 

• 

U 

loaloo  <a-tO'rt!f> 

1 

sOO«H 

fH 

X 

•  •  •  ••  •  •  • 

X 

»-4 

iC 

sotooiDO 

u 

>doac  <ooa^<o^ 

u 

H- 

•  •  •  #•  •  •  • 

c 

•J 

_l 

U 

X 

o 

u. 

4caj\0 

(J 

o 

O  C-  C 

o 

1  1  1  1  1  1  1 

o 

'2 

UitUtAj  lUUXftJUJ 

(D 

rr 

.  -  f*-  •'.  ’ 

tn  O'T  ^Si^CVjCM 

lU 

o 

n  mo  toxof^  o 

o 

Cl 

«  ••  ••••  • 

CM  CM*-*  CMf**  OJ  O 

o 

o 

z 

u 

fo  K)fO  rOKJrO  K)  fO 

UJ 

<M  CjCM  CV.CJCMCM  CM 

ui  ipinio  Si  -fi 

u 

UJ 

o 

X 

X 

►- 

t- 

►- 

CO 

4 

< 

u 

X 

u 

u 

>- 

4 

< 

u 

X 

®  (TNO  *^CJfO  ^iD 

u 

U 

a. 

X  Q<T' CTVTXT' 

Ut 

Ch'J'  ITKJ'CT'  J'  O' 

UJ 

UJ 

H* 

»  4--»  aa- 

X 

U. 

VJ 

88 


«  e« 

«  O  41 

«  • « 

41  *'41 

u. 

«  o« 

«  041 

XU. 

«  « 

41  « 

CO 

Ui 

«  « 

«  « 

V 

«  « 

4i  ♦ 

M3 

«  « 

«  « 

e 

«  « 

41  « 

<u. 

«  « 

CM 

XU. 

«  Ul** 

«  Ui4i 

o 

OUJ 

♦  -J« 

«  <-141 

«  u* 

«  U« 

CO 

«  >>« 

«  >-4l 

«  a« 

41  U.« 

' 

«  « 

«  41 

o 

«  « 

«  # 

M3 

«  <M« 

«  ^'-0 

in 

X 

♦  K)« 

«  o  41 

in 

«  <D« 

^- 

«  ^ « 

•ft  <4^-41 

CD 

«  « 

ft  ft 

• 

X 

«  « 

ft  '0 

< 

♦  a-« 

ft  CX.  ft 

II 

X 

«  Ui« 

ft  Ui.ft 

ID 

«  »-« 

ft  »-  ft 

U 

»-4 

«  co« 

ft  COft 

CO 

« o  « 

fto  ft 

♦  o  ♦ 

ftO  ft 

«  tf  1  « 

ftO  ft 

2 

u 

«  •  '  « 

ft  •  ft 

1 

« iO  ««« 

ft  CM  ^  ft 

O' 

X 

Ui« 

ftin  uift 

m, 

m  dK>  -9 

C/3 

X 

«  U4« 

ft  ^41 

X 

41  It 

ft  u.  >*ft 

A 

«o  « 

ftO  ft 

X 

«  x« 

ft  Xft 

u 

«o 

ft  CJ  ►•ft 

X 

u 

h- 

«  X  1  * 

ft  00  1  ft 

S 

c 

•J 

•  S  •«H4r 

ft  S  ft 

a: 

•J 

UJ 

♦  <o  c^  ♦ 

ft  C04-4  ft 

H- 

CD 

o 

«  XXI  « 

ft  X«4  ft 

O 

4>  •  % 

«  M 

(jJ 

u. 

X:  « 

ft  ft 

a. 

o 

«  gv  « 

41  CT  0«-4  ft 

«  #g>  « 

ft  •€>!  ft 

U3 

2 

^  *41 

ft  ft  ft 

o 

2> 

o 

*0^  u*  n 

ft  ft 

2 

•>4 

X 

«  3  X  « 

ft  O  Xft 

3 

X 

o 

«  o « 

ft  Oft 

w 

2 

Li 

«  »  iiky  « 

ft  pM  UI  ft 

a: 

«  ao  cj(^  « 

ft  oooco  ft 

o 

o 

«  •  4i 

ft  ft  ft 

1 

Ui 

41  ^  # 

ft  »  ft 

o 

••r-  « 

ft  N*  ft 

41  •  41 

ft  •  ft 

• 

w 

u 

4i  ^ 

ft  «  ft 

oe 

X 

*  KW'  ♦ 

ft  OO'  ft 

>- 

*  ^  Ck-4  « 

ft  JO*-«  ft 

< 

«  <  40  ■• 

ft  <o  ft 

CO 

41  7’  •  ♦ 

ftX  •  ft 

Ui 

«  o  in  ♦ 

ft  OCM  ft 

•J 

«  M  <9*  « 

ft  ^in  ft 

u 

•«  H-  X  « 

ft  X  ft 

u 

V 

4IU  MCJ« 

ft  o  II  cj  ft 

2 

o 

*  --Cl-O* 

ft  <^-o  ft 

CM 

X 

«  X 

ft  ft 

I 

o 

Q. 

«  U.  XCD4I 

ftULXOSft 

CD 

UJ 

41  M  41 

ft  •-*  ft 

1 

Ui 

41  .UZ  41 

ft»-.i2  ft 

Ul 

u. 

CO 

«  3C*-^4I 

ft  ft 

H* 

«  « 

ft  ft 

89 


f-B-2  61R-T0-&RCUN0  SPFCTRUM*  BASELINE 

T  hf  CRACK  at  the  BCGlANlMi  Of  BLOCK  20  C  =1,327377  02/06/81 
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AD-A'116  968  ROCKWELL  INTERNATIONAL  EL  SE8UND0  CA  NORTH  AMERICAN  —ETC  F/6  9/2 
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Appendix  A 


SUGGESTED  PROCEDURES  FOR  DETERMINING  CRACK -GROWTH -RATE  PARAMETERS 


To  use  the  CRKGRO  program  for  fatigue  crack  growth  life  predictions,  a 
set  of  crack -growth- rate  parameters  should  be  determined  for  each  material. 
The  following  paragraphs  briefly  describe  a  procedure  for  determining  those 
parameters . 


BASELINE  CRACK-GROVmi-RATE  CONSTANTS 

Fatigue  crack- growth- rate  constants  for  the  baseline  (constant-amplitude) 
crack-growth-rate  equations  used  in  CRKGRO  are  the  coefficient  C,  exponents 
n  and  m,  in  the  following  equation: 

da/dN  =  C  [(l-R)"*  K  for  R  >  0 

max 

plotting  the  dependent  variable  da/dN  against  the  independent  variable 
(l-R)"*  K__y  on  a  double  logarithmic  scale,  the  preceding  equation  represents 
a  straight  line;  i.e. , 

In  (da/dN)  =  In  C  +  n  In  [(1-R)'"  K^] 

where  n  is  the  tangent  of  the  angle  of  inclination  and  InC  is  the  distance 
of  intersection  with  the  In (da/dN)  axis  from  the  origin. 

The  procedures  to  deteimine  the  coefficient  C  and  exponents  n  and  m  of 
the  preceding  baseline  crack- growth- rate  equation  for  a  specific  material 
are  as  follows: 

1,  Perform  constant -amplitude  fatigue  crack  growth  test  per  ASW 
standard  E-647  or  equivalent  specification.  It  is  recommended  to 
perform  the  test  with  minimum  four  stress  ratios R=0,  0,3,  0.5  and 

0.7. 

2,  Convert  the  crack  size,  a,  versus  the  number  of  the  elapsed  cycles, 

N,  to  the  fatigue  crack  growth  rate,  da/dN,  using  any  of  the  two 
ASTM  recommended  data  reduction  techniques:  the  second  method  or 
the  incremental  polynomial  method. 

3,  Plot  the  independent  variable  [(l-R)""  K„^]  versus  the  dependent 
variable,  da/dN,  on  log-log  coordinates  with  assigned  m-values, 
either  by  hand  or  by  a  graphical  routine  such  as  PLOTRATE,  developed 
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by  Chang  et  al  (reference  11  ).  It  is  recommended,  in  general,  to 
input  the  following  set  of  m- values: 

m=0.4,  0.5,  0.6,  0.7,  0.8,  0.9  and  1.0 

Mote  that  if  m  =  1,  the  independent  variable  =  AK. 

4.  Examine  each  da/dN  versus  (l-R)*"  J^nax  plot-  Hio  m- value  of  a  plot 
which  shows  the  best  stress  ratio  layer  collapsing  is  the  proper 
value  to  be  used  in  the  analysis.  A  straight  line  is  then  drawn 
either  by  eye-balling  or  conputer,  with  least-square  routine  through 
all  the  collapsing  data  points.  The  intersection  of  this  straight 
line  with  the  ln(da/dN]  axis  gives  the  C-value,  while  the  tangent  of 
the  angle  of  inclination  of  the  straight  line  gives  the  value  of  n. 

A  typical  plot  generated  by  PLOTRATE  for  2219-T851  aluminum  is  shov.r 
in  figure  A-1.  It  shows  that  m  =0.6  has  the  best  collapsing  on  ai; 
the  data  points  at  various  stress  ratios  (R  =  0.01,  0.2,  0.3,  O.b 
and  0.7).  The  corresponding  C-  and  n-v^alues  are 

C  =  5.063  X  10  (in  ksi  unit) 

n  =  3.83 


LOAD  INTERACTION'  ?DDEL  PARAMETERS 


Fatigue  crack  growth  parameters  to  be  used  in  the  load  interaction 
model  built  into  CRKGRO  are  the  overload  shutoff  ratio,  Rg^,  the  acceleration 
index,  q,  and  the  threshold  intensity  factor  range  at  R  =  0, 


The  overload  shutoff  ratio,  Rgg,  and  the  threshold  stress  intensity 
factor,  ,  are  the  parameters  used  in  the  generalized  Willenborg  retard¬ 

ation  model  tor  calculating  the  retardation  coefficient,  <t'  ,  which  is 
defined  as;  / 

I  -Ik 


♦  = 


\ 


max 


thy 


<  ) 


R  -  1 
so 


where  Kjnax^h  threshold  stress- intensity- factor  value  corresponding  to 

the  value  in  the  following  form: 

‘-'max  ,  ■  -  W  -  »  ■  ■  « 

th  o 
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BASELINE  2219-T85:  Al'J)^  CRACK  GROWTH  RATE  C*."a  oWlkiToco 


ROCKWELL  EQUATION 


0  B-l-l  R  ■  0.010  X  9-1-2  R  •  0.200  g  9-2  R  .  0.300  Y  9-3  R  •  3  730 

♦  9-«-l  R  ■  0.010  •  9-6-2  R  -  0.010 

N  •  3.93  C  •  9.aB3E-IO  M  •  0.60 


Figure  Baseline  2219-T851  Aluminum  Crack-Clrowth-Rate  Data, 


To  determine  the  overload  shutoff  ratio,  a  series  of  single  overload 
variable-amplitude  cyclic  tests  are  required.  This  is  because,  for  a  given 
material,  the  overload  shutoff  ratio  is  dependent  upon  the  stress  ratio  and 
the  magnitude  of  the  compressive  stresses.  But,  to  be  practical,  it  is  rec¬ 
ommended  to  test  only  the  R  =  0  condition  to  obtain  the  Rgotoj  value.  This 
value  is  then  to  be  calibrated  using  a  typical  spectrum  which  has  similar 
characteristics  to  the  one  to  be  predicted. 

The  acceleration  index,  q,  is  the  exponent  in  Chang’s  crack-growth- rate 
equation  for  negative  stress  ratio,  which  is  expressed  as: 

da/dN  =  C  [(1  +  R^)^  K  ]”,  R  <  0 

max  * 

The  value  of  q  is  determined  using  the  following  relationship: 
q  =  [j{n(Y}/in(l  +  R^)]/n,  R  <  0 

where  y  is  the  ratio  of  the  crack  growth  rate  at  a  specific  negative  stress 
ratio  to  its  R  =  0  counterpart  obtained  from  test  data.  Hence,  based  on  the 
preceding  metholology  for  a  specific  negative  stress  ratio,  there  should  be  a 
specific  q  value.  However,  for  spectrum  loading  application,  it  is  not  very 
practical  to  generate  such  q-values.  The  average  q- value  approach  was  adopted 
by  CRKGRO.  The  average  q-value  for  a  material  can  be  selected  by  correla¬ 
tion  with  the  test  data  obtained  from  a  spectrum  loading  test  while  the  spec¬ 
trum  is  similar  to  the  one  to  be  predicted. 


OTHER  RARAMETERS  USED  IN  THE  METHODOLOGY 


There  are  several  other  parameters  needed  to  be  input  into  CRKGRO  for 
crack  growth  predictions,  including  the  cutoff  values  of  the  positive  and 
negative  stress  rations,  R^ut>  critical  values  of  the  stress  intensity 

factors  under  cyclic  loadings, 

R^ut  is  the  cutoff  value  of  the  positive  stress  ratio,  R^,  above  which 
the  material  is  assumed  to  have  no  further  stress  ratio  layering  effect  on 
the  crack  growth.  R^u^  is  the  cutoff  value  of  the  negative  stress  ratio,  R". 
below  which  the  material  is  assumed  to  have  no  further  acceleration  effect  on 
the  crack  growth.  These  values  are,  in  general,  able  to  be  determined  from 
multiple  stress  ratios  constant-amplitude  tests.  Yet,  for  the  spectrum  load 
application,  because  the  effective  stress  ratio  is  used  in  the  load  interac¬ 
tion  model,  again,  the  calibrated  value  should  be  obtained  using  the  spectrum 
test  data. 
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There  are  two  types  of  Kcr  values:  l^cr(a)  ^cr (g) '  *^cr(a) 

critical  value  of  the  stress  intensity  factor  under  cyclic  loading  at  the 
maxinun  depth  point  A  of  a  PTC.  In  general,  the  value  of  K^-i-(a)  considered 
approximately  equal  to  material  plane-strain  toughness.  K^-rCc)  critical 

value  of  the  stress  intensity  factor  under  cyclic  loading  at  the  maximum  length 
point  c  of  a  PTC  or  TC.  TTie  value  of  considered  approximately  equal 

to  material  plane-stress  toughness. 
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Appendix  B 

CRKGRO  STRESS  ri^TENSITY  FAC1UKS  LIBRARY 


4  ^ 
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CUKGKO  Stress  Intensity  factor  Library 


CRKGRO  Stress  Intensity  Factor  Library 


CRKGRO  Stress  Intensity  Factor  Library 


no 


where  Q  =  1  +  1.464  (c/a) 


CRKGRO  Stress  Intensity  I'actor 


CRKGRD  Stress  Intensity  Factor  Library 


wfiere  Q  -  1  +  1.464  (c/a) 


ess 


CRKGRO  Stress  Intensity  Factor  Library 


where  Q  =  1  *•  1.464  (c/n) 


•actor 


ty  lac tor 


t  ivc  Tfst  Mctlioil  for  (!onsi 
78  Amuuil  Bo<,)k  of  AS'l'M  St; 


Q^KCRO  Stress  Intensity  Factor  Library 


*Reference:  Nevvinan,  ,J.  C.,  "Predicting  [ailnrc  of  Specimens  with  liithei'  Surface  Cracks  or  Cornci'  Cracks  at 
Holes,"  NASA  TN  I)  8241,  June  p  7 


CRKGRL')  Stress  Intensity  I-actor  I  ibrary 


''I’rediet  ing  lailiire  of  Specimens  with  iiithcr  Surface  (backs  or  Corner  (  rac 


Kci crenc 


Appendix  C 

TYPICAL  FIGKTER  AIR-TO-GROUND  BASELINE  MISSION 
SPECTRUM  TABLE 
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Appendix  C 


EXPERIMENTAL  VERIFICATION  PROGRAM,  TEST  GROUP  I -A 

FIGHTER  BASELINE  TEST  ^ 

TEST  F-B-2,  AIR-TO-GRDUNB  MISSION,  DESIGN  LIMIT  STRESS  (DLS)  »  30  KSI 
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EXPERIMEVT-AL  VERIFICATION  PROGRAM,  TEST  GROUP  I-A 
FIGHTER  BASELINE  TEST  (CONTINUED) 

TEST  F-B-2,  AIR-TO-GROUND  MISSION,  DESIGN  LIMIT  STRESS  (DLS)  =  30  KSI 
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FIGHTER  BASELINE  TEST  (CONTINUED) 

TEST  F-B-2,  AIR-TO-GROUND  MISSION,  DESIGN  LIMIT  STRESS  (DLS)  -  30  KST 
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